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ABSTRACT 

Detection of the signature of Rayleigh scattering in the transmission spectrum of 
an exoplanet is increasingly becoming the target of observational campaigns because 
the spectral slope of the Rayleigh continuum enables one to determine the scaleheight 
of its atmosphere in the absence of hazes. However, this is only true when one ignores 
the refractive effects of the exoplanet’s atmosphere. I illustrate with a suite of simple 
isothermal clear Jovian H 2 -He atmosphere models with various abundances of water 
that refraction can decrease significantly the spectral slope of the Rayleigh continuum 
and that it becomes flat in the infrared. This mimics a surface, or an optically thick 
cloud deck, at much smaller pressures than one can probe in the non-refractive case. 
The relative impact of refraction on an exoplanet’s transmission spectrum decreases 
with atmospheric temperatures and increases with stellar temperature. Refraction is 
quite important from a retrieval’s perspective for Jovian-like planets even at the high¬ 
est atmospheric temperatures (1200 K) considered in this paper, and for all stellar 
spectral types. Indeed, refraction breaks in large part the retrieval degeneracy be¬ 
tween abundances of chemical species and the planet’s radius because the size of 
spectral features increases significantly with abundances, in stark contrast with the 
non-refractive case which simply shifts them to a larger or smaller effective radius. 
Abundances inferred assuming the atmosphere is cloud-free are lower limits. These re¬ 
sults show how important it is to include refraction in retrieval algorithms to interpret 
transmission spectra of gas giants accurately. 

Key words: atmospheric effects - methods: numerical - planets and satellites: at¬ 
mospheres - planets and satellites: gaseous planets - radiative transfer. 


1 INTRODUCTION 


Transmission spectroscopy is a powerful observational 
method that allows one to probe the atmosphere of a tran¬ 
siting exoplanet. Occultation by the body of the planet and 
absorption by the planetary atmosphere cause a drop in 
the stellar flux during the planetary transit. The absorption 
cross section of the exoplanet as a functi on of wavele ngth, 
expressed in terms of its effective radius llBrow nl l200lll , de¬ 
pends on the transmission of the atmosphere at the termi¬ 
nator as a function of altitude, as well as the radius of the 
planet at a reference atmospheric pressure. 

To understand which atmospheric parameters can 
be retrieved from the ana l ysis o f transmission spectra, 
Ibecavelier des Etangs et al.l ( 20081 ') derived a few simple 
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relationships for an isothermal atmosphere with a uni¬ 
form vertical composition. They showed numerically that 
the effective radius of an exoplanet occurs at an alti¬ 
tude where the integrated optical depth along the line of 
sight grazing the planetary limb, or slant optical depth, 
is about 0.56. This was la ter demonstrated analytically 
by de Wit fc Seager Usi ng this simple relation, 

Ibecavelier des Etangs et al.l (l2008l l also showed that the size 
of spectral features scales with the atmospheric scaleheight. 
Assuming H 2 Rayleigh scattering to be the most impor¬ 
tant source of extinction from 0.6 to 1 p.m, they derived 
the scaleheight, and hence the temperature, of HD 189733b 
from transit data obtained with the Hubble Space Telescope 
JPont et al.ll200^ '). 


Since the spectral dependence of Rayleigh scattering 
cross sections is well known, spectral regions dominated 
by Rayleigh scattering are particularly useful to deter- 
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mine the scaleheight of an atmosphere by measuring the 
change in the effective planetary radius with wavelength 
dLecavelier des Etangs et'aI]l20oA iBenneke fc Seagejll2012l : 
iHowe fc Burrowd 20121 ). This Rayleigh slope can be used 
to discriminate between various bulk atmospheric com¬ 
position, provided the atmosphere is relatively clear of 
hazes or clouds. In the presence of extended vertical hazes, 
the transmission spectrum will instead bear the signa¬ 
ture of the hazes’ optical properties (iHow e fc Bu rrowsI 
I2OI2I : iRobinson et al.l I2OI4I : fwakeford fc Sin'S " 2015l l. Ei¬ 
ther way, this slope provides an important constraint 
on the properties of exoplanetary atmospheres, and has 

been detected in _ a f ew exoplan e ts, e .g. HD 209458b 

I Knutson et al. 200?! : ISing et ^ 2008 ll. HD 189733 b 


I Pont et al.l 200^, HA T-P-5b | Southwo rth et al] |2012|), 


WAS P-12b llSing et a l.l 2013ll. GJ 3470b dNasc imbeni et al.l 
l2013h . WA SP-6b djorda n et al. [20131 : iMkolov et ^ 20151)^ 
QatM -2b (iMancini et aL~| 2014|), WASP- 31b ( Sing et ^ 
l2015h . and WASP-103b I Southworth et alll2015ll . 


In spectral region where molecular absorption domi¬ 
nates, the atmospheric scaleheight can still be estimated 
from the size of absorption features if one assumes that 
molecules are well-mixed in the atmosphere. However, it 
is very difficult to determine the absolute abundances of 
species and the pressure regions that are probed, as degen¬ 
erate solutions exist which can match observations. Indeed, 
changes of several orders of magnitude in molecular abun¬ 
dances can be compensated with a small change in the plan¬ 
etary radius at the reference pressure (iGriffithI |2014^ . Al¬ 
though the absolute mole fraction of species is ambiguous, 
their relative abundances are well const rained through the 
relati ve strength of absorption features (IBenneke fc SeagCT 
l 2012 li in a haze-free atmosphere. 

However, none of these previous investigations of the 
retrieval of atmospheric parameters considers the refrac¬ 
tive effects of the exoplanetary atmosphere on its trans¬ 
mission spectrum, in stark contrast with stellar occulta- 
tions. The use of refractive effects to probe planetary atmo¬ 
spheres has a long standing history in Solar system explo¬ 
ration which started wi th the occultation of a Arietis by the 
atmosphere of Jupiter (iBaum fc Codel [19531 1. This pioneer¬ 
ing work paved the way for stellar occultation by planetary 
atmospheres to become a powerful method by which one 
can derive not only the density vertical profile of planetary 
atmospheres, but also deduce the vertical profile of minor 
chemica l species which absorb stellar radiation (see the re¬ 
view bv ISmith fc Huntenlll99(il '). However, the observational 
geometry of a stellar occultation is different from that of a 
transiting exoplanet. Indeed, in a stellar occultation, the ra¬ 
diation source is infinitely far away compared to the planet 
and the observer, whereas in an exoplanet transit geometry 
the observer is infinitely far away from the exoplanet and its 
star. 

iHubbard et al.l 1 I 2 OOIII first investigated the impact 
of refraction on exoplanet transmission spectra and de¬ 
duced that refractive effects were negligib l e for the hot 
Jupiter HD 209458b. However. iHui &: Seagerl (l2002f) adapted 
gravitational microlensing theory to atmospheric refrac¬ 
tion and showed that refracti on could b e imp ortant 
for some transiting exoplanets. ISidis fc Saril (l 2010 l 'l used 
a more conventional ray-tracing approach and came to 
the same conclusion. They both found that even in 


the absence of absorption, an exoplanet’s effective ra¬ 
dius is still quite large because the higher density re¬ 
gions of planetary atmospheres deflect stellar radiation away 
from the observer. Absorption by chemical species were 
later combined with ray-tracing JCarcfa Munoz et al.ll2M2 : 

iBetremieux fc Kalteneggerl I2OI3I : Misra, Meadows & Crisp 

2014) to compute transmission spectra of an Earth-Sun ana¬ 
log, and showed that refraction decreases signihcantly the 
contrast of H 2 O features because transmission spectroscopy 
cannot probe most of the troposphere, which holds the bulk 
of E arth’s water inventory. _ 

iBetremieux fc Kaltenegg^ ll2014l ') explored the param¬ 
eter space of an Earth-like planet receiving the same stellar 
flux as Venus, Earth, and Mars, for different spectral type 
host stars along the Main Sequence. They showed that one 
can only probe an Earth-like planet with the same insolation 
as Earth to a pressure of 1 bar if the planet orbits a M5 or 
cooler star, assuming it maintains a similar temperature- 
pressure profile. As the temperature of the host star in¬ 
creases, the contrast of spectral features decreases, poten¬ 
tially resulting in flat spectral regions (see their Figure 11). 
Refraction also produces wings brighter than the contin¬ 
uum just outside of transit dHui fc Seaget|[200^ : ISidis fc Saril 
l2010l : [^rcfa Munoz et al.l|2012ll because some of the stellar 

radiation is deflected toward the observer while the star is 
unocculted. These brighter wings are sensitive to th e aerosol 
content of an atm osphere ([Misra fc Meadowj 201^). 

More recently. IBetremieux fc Kaltenegger ( 2015lf devel¬ 
oped new analytical expressions to describe the integrated 
column abundance along the curved path of a refracted ray, 
as well as the total deflection of the ray. They discovered 
that thick atmospheres form a thin refractive boundary layer 
where the column abundance and its deflection, go from fi¬ 
nite to infinite values as the grazing altitude of a ray de¬ 
creases and approaches a lower refractive boundary. Indeed, 
at this lower boundary, the radius of curvature of the ray in¬ 
duced by refraction matches the radial position at which the 
ray grazes the atmosphere, and the ray follows a closed cir¬ 
cular path. These results diverge substantially from expres¬ 
sions developed for stellar occultations (Baum & Code 1953) 
and still in use today. Rays that penetrate below this lower 
boundary will spiral into the atmosphere and be absorbed. 
Hence, atmospheric regions located below this lower bound¬ 
ary cannot be probed by techniques based on the transmis¬ 
sion of radiation through an atmosphere, irrespective of the 
planet-star geometry. Furthermore, the effective scaleheight 
of the atmosphere, i.e. as perceived by an observer, decreases 
dramatically as this lower boundary is approached. 

What impact does refraction, and in particular this new 
theory, has on transmission spectroscopy of giant exoplan¬ 
ets? Does it affect the Rayleigh scattering slope, or the re¬ 
trieval of atmospheric parameters? This paper explores these 
issues for Jupiter-sized planets with a Jovian composition 
across various planetary temperatures and spectral type of 
the host star. 


2 DETAILS OF SIMULATIONS 

To illustrate the impact of refraction on transmission spec¬ 
troscopy of giant planets and on the retrieval of their 
atmospheric properties, I compute a suite of transmis- 
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sion spectra from 0.4 to 5.0 [rm for a Jovian-like planet 
with a H 2 -He atmosphere orbiting stars of different spec¬ 
tral type along the Main Sequence (M2, K5, G2, FO). I 
consider isothermal spherically-symmetric well-mixed at¬ 
mospheres for three different temperatures (300, 600, and 
1200 K), and various constant mole fractions of water 


vapour (0, 10 


10 " 


10 " 


10 " 


and 10 ). These 


simulations are done both with and without refraction. 
Stellar limb-darkening is ignored in these simulations, just 
as it was ignored in all derivations of the Rayleigh slope 
llLecavel ier des Etangs et al. 20081: ^nneke fc Seagerll2012l : 
iHowe &: BurrowsI I2OI2I : Ide Wit fc Seageil 120131 '). to have a 
common reference frame and enable a direct comparison 
with previons simulations of transit spectra. Mnch of the 
theory and methods at the hear t of these simulations have 
alrea dy been described in detail (iBetremieux fc Kalteneggerl 
l2013l . 2014, 2015), bnt are summarised in this section for the 
convenience of the readers. 


2.1 Planetary atmospheres 

I bnild-np the density vertical profile of an isothermal home- 
geneous atmospheres in hydrostatic equilibrium with a grav¬ 
itational acceleration which varies with altitnde, for three 
different temperatures (300, 600 , and 1200 K), in the man¬ 
ner de scribed in Section 3.3 in IBetremieux fc Kaltenegg^ 
ll2015h . The radius {Rp = 69911 km) and mass (Mp = 
1.8986 X 10^^ kg) of the exoplanet are the same as 
Jupiter’s. The atmosphere has a He mole fraction of 0.1357, 
as measured on Jup i ter b y the Galileo probe (listed in 
iLodders fc Feelev Jr.l Il998h . while the remainder of the 
atmosphere is predominantly composed of H 2 and trace 
amounts of water. The planetary radius (Rp) is dehned at 
a reference pressure of 1 atm (atmosphere), and the highest 
pressure considered in the simulations, or “surface” pressure, 
is 500 atm. The densities are computed on a fine computa¬ 
tional grid with thickness Azcomp, and then sampled on a 
coarser grid with altitude intervals Azaamp- Table [T] lists 
the thickness of the layers making the computational and 
sampling grids, along with the chosen atmospheric thickness 
{AZatm), and the resulting “surface” altitude (za) for each 
temperature. The atmospheric thickness defines the radius 
of the top atmospheric boundary (Rtop), given by 

Rtop ~ Rp Za AZatm,’ (1) 


2.2 Planet-star distance 


To determine the distance between the planet and its stel¬ 
lar host, I assume that the isothermal temperature of the 
planet (Tp) matches the mean planetary emission tempera¬ 
ture and that the planet has a circular orbit. Under these 
assumptions, the semi-major axis (a) of the planet is given 
by, 


VI — As 



( 2 ) 


where Ap is the Bond albedo of the planet, while T* and 
7?* are the stellar effective temperature and radius, respec¬ 
tively. Table[2]lists these stellar parameters for the four stel¬ 
lar spectral type (M2, K5, G2, FO) along the Main Sequence 


Table 1. Input and derived model atmosphere 
parameters for different atmospheric tempera¬ 
tures 


Input 

Tp (K) 

300 

600 

1200 

(km) 

1370 

2760 

5640 

^'^samp (km) 

5 

10 

20 

^'^comp (km) 

0.01 

0.02 

0.04 

Derived 

Za (km) 

-260.70 

-519.48 

-1031.28 


Table 2. Stellar parameters llCoxll200(ili 


Spectral type 

T* (K) 

R* ( 72 ©) 

FO 

7300 

1.50 

G2 

5778^ 

1.00 

K5 

4410 

0.72 

M2^ 

3400 

0.44 





“ iLodders fc Feglev Jr.l lll998l ') 

^ values from iReid fc: Hawlevi ll20n.5l'l 


considered in this paper. A Bond albedo of 0.30 is chosen, 
similar to those of the giant planets in our Solar system. 


2.3 Column abundances and ray deflection 

To compute the transmission spectrum of an exoplanet, one 
must determine the altitude-dependent incremental column 
abundance encountered by a ray along its curved refracted 
path, as well as the total deflection of the ray, as a function 
of a ray’s grazing height (ro), i.e. planetary radius + grazing 
altitude of the ray. I comp ute these quantities with MAKEX - 
OSHELL (see Section 3 in iBetremieux fc Kalteneggerll2015l l. 
including the contribution from atmospheric regions above 
the top atmospheric boundary, with 80 rays spread uni¬ 
formly in grazing altitude between the top and bottom 
boundaries. These grazing altitudes also define the bound¬ 
aries of the atmospheric layers used in the computation. For 
the refractive case, I use a refractivity at standard temper¬ 
ature and pressure (vstp) of 1.23x10”^, that of the Jovian 
atmosphere at 1 qm, and a refractivity of 0 otherwise. The 
location of the bottom atmospheric boundary depends on 
the condition of the simulations. In the refractive case, the 
bottom boundary can be either the critical or the lower re¬ 
fractive boundary (see Section 12.411 . otherwise the bottom 
boundary is the “surface”. 


2.4 Critical and lower refractive boundaries 

The lower refractive boundary is located at a height 
(r) where, in a spherically symmetric atmosphere, a 
grazing ray follows a circular path (see Section 2.2 in 
IBetremieux fc Kalteneggerl I 2 OI 5 I I. Atmospheric regions be¬ 
low this boundary cannot be probed by exoplanet trans¬ 
mission spectroscopy or by stellar occultations because rays 
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Table 3. Refractive boundaries of Jupiter-sized planets with different isothermal tempera¬ 
tures orbiting stars of various spectral type (see Section l2.4|l 


T (K) 

Critical Boundary 

300 

n (amagat^ 

P (bar) 

600 

n (amagat) 

P (bar) 

1200 

n (amagat) 

P (bar) 

FO star 

0.33 

0.37 

1.78 

3.97 

8.64 

38.45 

G2 star 

0.54 

0.60 

2.78 

6.20 

12.21 

54.34 

K5 star 

0.92 

1.02 

4.45 

9.90 

16.25 

72.35 

M2 star 

1.56 

1.73 

6.62 

14.74 

18.75 

83.45 

Lower Boundary ^ 

4.88 

5.43 

9.74 

21.69 

19.38 

86.27 


^ 1 amagat = density at standard temperature and pressure (STP) 
^ independent of stellar spectral type 


that reach these regions spiral deeper into the atmosphere 
until they are absorbed. This boundary is located at a height 
where 


(_j4r)_ \ ^ 
\l + v{r)) H 


(3) 


and this depends only on the properties of the atmosphere. 
Here, v{r) is the height-dependent refractivity of the atmo¬ 
sphere, and H is the density scale height. 

In reality, one can only p robe to the critical refrac¬ 
tive b oundary (see Section 2.2 in [Betremieux fc Kalteneggeij 
I 2 OI 4 I '). This boundary occurs at a density where the atmo¬ 
sphere deflects a ray coming from the opposite stellar limb 
toward the observer. Hence, this boundary depends both on 
the properties of the atmosphere and on the geometry of the 
planet-star system, i.e. on the angular size of the star as seen 
from the planet. When the planet occults the central region 
of its host star, this boundary is located at the same density 
all around the planetary limb, where rays are deflected by a 
critical deflection (tJc). During the course of a transit, this 
symmetry is broken and one can probe higher densities on 
one side of the planet, and lower densities on the opposite 
side. However, since the ray deflection tends to infinity as 
the grazing height of the ray approaches the lower boundary, 
the critical boundary can only asymptotically approach but 
never reach the lower boundary. 

MAKEXOSHELL computes the quantity on the left- 
hand side of equation [3] on the fine altitude computational 
grid, and determines the lower boundary by choosing the 
altitude for which this value is closest to but less than one. I 
also use MAKEXOSHELL to find the location of the critical 
boundary by computing the ray deflection for a grid of graz¬ 
ing heights, and then determine which one has the smallest 
difference with the critical deflection of the planet-star sys¬ 
tem, which is given by 


sinwc 


ro -I- R* 
a 


(4) 


Table [3] shows the critical boundaries for the four spectral 
type stars from Table [2] both in terms of the density and 
pressure of the atmosphere, for the three isothermal temper¬ 
atures considered. It also lists the lower boundaries which 
are independent of the spectral type of the host star. 


2.5 Radiative transfer 

The radiative transfer code jBetremieux fc Kalteneggerl 
I 2 OI 3 II combines the MAKEXOSHELL-computed column 
abundances {Nik) of the atmospheric layer encountered 
by the ray, with the average mole fraction (fijk) and 
extinction cross section {aj) of the chemical species, to 
compute the integrated optical depth (n) along the ray, 
using 

ri = '^aikNik, (5) 

k 

where 


^ik — ^ ^ ^j fijk (h) 

3 

is the average extinction cross section of the layer along 
the ray. Since all the molecular species are well-mixed 
in these simulations, fijk is independent of the atmospheric 
layer and of the ray, and is simply the average mole fraction 
of the species. The transmission (Ti) of the ray is 
simply given by 


Ti = e- 


(7) 


Earlier versions of the code (see e.g. I Johnson et al.l Il995l : 
iKaltenegger fc TraubI l2009l f added thermal emission to 
the transmitted stellar radiation. However, for a one¬ 
dimensional spherically symmetric atmosphere, the thermal 
emission contribution to an exoplanet light curve is constant 
throughout the transit, and the calculated transit depth is 
indep endent of thermal emission . To sa ve on computation 
time, iBetremieux fc Kaltenegg^ l|201,3h modihed the code 
so that one can choose not to include thermal emission. 

Although the line opacity d atabase is identical to that 


in Kaltenegger fc TraubI (l2009l L IBetremieux fc KalteneggeH 
1 I 2 OI 3 I 1 replaced the original Rayleigh scattering parametri- 
sation, appropriate only for Earth, with a database of 
Rayleigh scattering cross sections for N 2 , O 2 , CO 2 , and Ar 
(see their Table 2). They also included new routines to com¬ 
pute optical depths from continuous absorption cross sec¬ 
tions, and a new database of ultraviolet and visible absorbers 
(see their Table 1). 

For this paper, I have added the Rayleigh scatter¬ 
ing cross section of molecular hydr ogen and helium to the 
database. I use the formulation by iFord fc Brownel (Il973fl 
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for the H 2 cross section as a function of its rotational quan¬ 
tum number J, from which I compute the cross section for an 
equilibrium H 2 population at 300 K. Since the temperature 
dependence of the cross section is very small (e.g. about 0.6 
per cent change between 273.15 and 2000 K at 1 pm), I use 
the 300 K cross section for all temperatures. For helium, 
I compute its STP refr activity with the e xpression from 
iMansfield fc Pecl3 (Il969l ~) , also described in IWeberj ll2003ll , 
and obtain its Rayleigh cross section {(Ju) with Rayleigh’s 
formula: 




327r® 

3 


\ 2 

VSTP \ 4 

- W Fk- 

nsTP ) 


(8) 


Here, nsTP is the number density at standard temperature 
and pressure, also known as Loschmidt’s number, w is the 
wavenumber of the radiation, and Fk is the King correction 
factor. For helium, Fk is equal to one at all wavenumbers. 

The only other source of opacity considered in these 
simulations is that of water vapour. The algorithm for the 
line-by-line computa tion of optical depths is described in 
iTraub fc Stietl (Il976l l. To reduce the computation time, only 
H 2 O lines which have a peak optical depth greater than 
10“^ within an atmospheric layer are included, and the line 
profile is computed to an optical depth of 10“® or to the 
edge of the spectral region, whichever is closer to the centre 
of the line. The atmospheric transmission is computed from 
1800 to 26000 cm“^ every 0.05 cm“^, and then averaged 
in 4 cm“^-wide bins. I choose a spectral region which is 
larger than that of interest to mitigate edge effects. The 
transmission spectrum is eventually displayed with respect 
to wavelength from 0.4 to 5.0 qm. 

To compute the effective radius (Re//) of the transiting 
exoplanet, I use the 80 rays from MAKEXOSHELL (see Sec¬ 
tion [Oil to partition the atmosphere into 80 annuli evenly 
spread in altitude. I then map these altitude boundaries to 
their corresponding impact parameter {bi) across the plane¬ 
tary limb, which is what an observer perceives. The impact 
parameter is given by 


= (1 + iy{roi))roi. 


(9) 


Taking the transmission of each annulus to be the average 
of the transmission at the boundaries, and ignoring stellar 
limb-darkening, the effective radius is then given by 



(|Betremieux fc Kalteneggeil I 2 OI 3 I 1 where N is the number 
of annuli. The {N + ray is that which grazes the top 
atmospheric boundary, chosen so that both atmospheric ab¬ 
sorption and deflection are negligible. Thus Tjv-i-i = 1, and 
6jv-i-i = Rtop- Atmospheric regions with densities higher 
than the density at the lower atmospheric boundary (i.e. 
the “surface”, lower refractive boundary, or critical refrac¬ 
tive boundary depending on the simulation) are implicitly 
assumed to be opaque at all wavelengths. One advantage to 
including refraction is that one does not need to compute 
the transmission of atmospheric regions located at higher 
densities and pressures, or lower altitude, than the lower re¬ 
fractive boundary. I display the transmission spectrum in 
terms of the effective atmospheric thickness (Az^ff) of the 
transiting exoplanet with respect to the 500 atm “surface”, 


which I compute with 

~ Re// — (Rp + Zs). (11) 


3 RESULTS AND DISCUSSION 


3.1 Rayleigh scattering slope 

To explore how refraction affects the Rayleigh scattering 
slope, I first compute transmission spectra where the only 
source of opacity is Rayleigh scattering from H 2 and He. 
Figures [T] through [3] show the 0.4-5.0 gm H 2 -He Rayleigh 
continuum under various conditions for the 300, 600, and 
1200 K isothermal atmospheres, respectively. The legends 
in all three figures are identical. The 500 atm “surface” is 
located at 0 km altitude, and the 1 atm pressure level is 
indicated for reference. Since the main effect of temperature 
is to change the scale height of the atmosphere, I scale the 
range of values on the ordinate with the planetary tempera¬ 
ture so that approximately the same number of scaleheight 
are displayed in all figures. It is not exactly the same be¬ 
cause the gravitational acceleration varies with altitude in 
the simulations, and so does the scaleheight. It is useful to 
plot things this way to notice secondary effects subtler than 
the primary effect of the scaling of spectral features with 
scaleheight. Indeed, even if transmission spectra have the 
same aspect from one figure to the next, the actual size of 
spectral features still differ by roughly the ratio of the plane¬ 
tary temperatures. I have done this in all figures throughout 
the paper. 

Without refraction, the transmission spectrum of a 
Rayleigh scattering isothermal atmosphere should fol¬ 
low a straight line when plotted in terms of the ef¬ 
fective atmospheric thi ckness versus the natural loga¬ 
rithm of the wavelength ( Lec avelier des Etangs e t al.ll2008l : 
iBenneke &i SeageH I 2 OI 2 I : Howe fc Burrows 2012 ), as long 
as the Rayleigh scattering cross section follows a constant 
power law with wavelength. The Rayleigh continuum for 
the non-refractive case (shown by the short-dashed line in 
all three figures) has roughly a constant slope except long- 
ward of 3 p.m where the Rayleigh cross section is suffi¬ 
ciently weak that the opaque regions below the “surface” 
contribute significantly to the transmission spectrum. In 
this case, these regions introduce a significant step function 
in the altitude-dependent atmospheric transmission profile, 
thus violating the conditions und er which the result from 
iLecavelier des Etangs et al.l ll2008li holds, and changing the 
Rayleigh slope. 

A noticeable secondary effect is that the spectrum also 
shifts to lower altitudes with increasing temperatures. This 
occurs because I refer ence the planetary radius to a pressure, 
as is often done (e.g. IBenneke fc Seagem2012l : Swain, Line, 
& Deroo 2014), rather than to a density. As temperature 
increases, the same pressure is achieved with lower densi¬ 
ties. Thus the densities required to achieve a slant optical 
thickness of 0.56 are located at lower altitudes and the trans¬ 
mission spectrum has an extra downward shift, beyond the 
primary scaleheight-induced factor, as the planetary temper¬ 
ature increases. It would seem that referencing the planetary 
radius to a density would be a better choice as the transmis¬ 
sion spectrum would not shift up or down with temperature. 
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0.4 0.7 1 2 3 4 5 

VACUUM WAVELENGTH (/lim) 

Figure 1. H 2 -He Rayleigh scattering transmission spectra of a 
Jupiter-sized exoplanet with a 300 K isothermal Jovian atmo¬ 
sphere, expressed in term of the atmosphere’s effective thick¬ 
ness above the 500 atm pressure level (as for all figures). The 
short-dashed line shows the spectrum without refraction, while 
the coloured triple-dot-dashed lines show spectra with refraction 
for different Main Sequence spectral type of the host star. The 
long-dashed line shows the spectrum with refraction if one could 
somehow probe the atmosphere to the lower refractive boundary, 
irrespective of the planet-star distance. The thick horizontal solid 
line shows the altitude of the 1 atm pressure level, for reference 
(shown in all figures). See Table |3] for the density and pressure 
levels of the critical and lower refractive boundaries. 



0.4 0.7 1 2 3 4 5 
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Figure 2. H 2 -He Rayleigh scattering transmission spectra of a 
Jupiter-sized exoplanet with a 600 K isothermal jovian atmo¬ 
sphere. See caption in Fig. [T] for more details. 

but I nevertheless adopt this common convention in this pa¬ 
per. 

In the refractive case, this simple picture changes dra¬ 
matically. Indeed, the black long-dashed line shows the 
transmission spectrum if one can somehow probe to the 
lower refractive boundary, irrespective of the planet-star dis¬ 
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Figure 3. H 2 -He Rayleigh scattering transmission spectra of a 
Jupiter-sized exoplanet with a 1200 K isothermal jovian atmo¬ 
sphere. See caption in Fig. [T] for more details. 

tance, which produces the smallest differences between the 
refractive and non-refractive cases. The infrared transmis¬ 
sion spectrum is flat, similar to that of an optically thick 
cloud or a surface, and gradually merges with the non- 
refractive case at lower wavelengths. The wavelength above 
which the Rayleigh slope deviates significantly from the non- 
refractive case increases with temperature, and changes from 
about 0.6 pm at 300 K to about 1 pm at 1200 K. The frac¬ 
tional change of the spectrum with and without refraction 
increases with wavelength and decreases with temperature, 
and it is still quite important even at 1200 K above 3 pm. 

Accounting for the location of the critical boundary as¬ 
sociated with the spectral type of the host star increases 
the effective thickness of the atmosphere and changes the 
Rayleigh slope further (see triple-dotted dashed coloured 
lines in Figures [T] through [3)). This correction can be quite 
important at the cooler planetary temperatures (see Fig¬ 
ure [TJ. As the temperature of the planet increases, the dif¬ 
ferences between the spectra associated with a FO and a 
M2 star decrease. At 1200 K (see Figure |3]), there is hardly 
any differences between the transmission spectrum around 
a FO and a M2 star, and both are very similar to the spec¬ 
trum where one can probe to the lower refractive boundary. 
Indeed, as the planetary temperature increases, the critical 
boundary approaches the lower refractive boundary. Since 
the effective atmospheric scaleheight decreases dramatically 
at the appro ach of the lower refractive bound ary (see Sec¬ 
tion 4.2.3 in [Betremieux fc Kalteneggerj[2015h . the altitude 
separation between the different critical boundaries also de¬ 
creases. 

Hence, the discovery of the re fractive boundary layer by 
Betremieux fc Kaltenegg^ (l2015l ) turns out to be extremely 
important for the transmission spectra of giant exoplanets, 
because even in the absence of clouds or hazes, one cannot 
probe the atmosphere to arbitrarily large pressures in spec¬ 
tral regions between molecular bands. Indeed, this refractive 
boundary layer not only mimics a surface but also decreases 
substantially the effective scaleheight of the atmosphere near 
that boundary. This, in turns, modifies not only the effec- 
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Figure 4. Transmission spectra without refraction of a Jupiter¬ 
sized exoplanet with a 300 K isothermal Jovian atmosphere for 
various vertical homogeneous mole fractions of water (see inset 
legend). The water-free Rayleigh scattering transmission spectra 
from Fig. [J are overlaid for comparison (see caption and inset 
legend in Fig. [T] for more details). Refraction will suppress water 
features which lie below the refractive Rayleigh curve associated 
with the host star’s spectral type (coloured triple-dotted dashed 
lines). 


Figure 5. Transmission spectra with refraction of a Jupiter-sized 
exoplanet, with a 300 K isothermal Jovian atmosphere for vari¬ 
ous vertical homogeneous mole fractions of water (see inset leg¬ 
end), if one could somehow probe the atmosphere to the lower 
refractive boundary, irrespective of the planet-star distance. The 
non-refractive Rayleigh continuum (short-dashed curve) is shown 
for comparison. 


Live radius of the Rayleigh continuum of the hotter exoplan¬ 
ets, but also decreases the Rayleigh scattering slope of the 
transmission spectrum. If one wants to infer th e atmospheric 
scaleheight (iLecaveli er dej Etangs et^_ahl l200^ or the mass 
of the exoplanet ( de Wit fc SeageJ 20131 ) from transmission 
spectra, one must be careful to evaluate the Rayleigh scat¬ 
tering slope in spectral regions that are largely unaffected by 
refraction, given the planetary temperatures. However, even 
then, the Rayleigh slope might still be affected by clouds or 
hazes, depending on their vertical distribution. 


3.2 Abundance retrieval 

To explore how refraction affects the retrieval of atmospheric 
abundances of chemical species, I now add various mole frac¬ 
tions of vertically-mixed water vapour (10“^^, 10“^°, 10~®, 
10“®, and 10”“^) to the simulations. I use water vapour as an 
example because its spectral signature is often sought after 
in exoplanetary transmission spectra, but the discussion in 
this section is generally applicable to other molecules (e.g. 
CH 4 , NH 3 , CO, CO 2 ) as well. 

Figure 0] shows the resulting transmission spectra 
(coloured solid lines) in the non-refractive case for the 300 K 
isothermal atmosphere compared with the various Rayleigh 
continua from Figure [T] Note that a H 2 O mole fraction of 
10“^^ does not produce any noticeable absorption features 
above the Rayleigh continuum and is not shown. Changing 
the water abundance mostly shifts features up or down with 
respect to the non-refractive Rayleigh continuum (short- 
dashed line). In spectral region where molecular absorp¬ 
tion is not much stronger than Rayleigh scattering, i.e. 
from 0.7 to 2 pm for H 2 O, the size of absorption features 
are pretty sensitive to the abundance of the responsible 
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Figure 6. Transmission spectra with refraction of a Jupiter-sized 
exoplanet, with a 300 K isothermal Jovian atmosphere for various 
vertical homogeneous mole fractions of water (see inset legend), 
orbiting a FO star. The non-refractive Rayleigh continuum (short- 
dashed curve) is shown for comparison. 


species. However, with the rapid decrease with wavelength 
of the Rayleigh scattering cross section, water abundances 
can change by several orders of magnitude above a certain 
threshold with hardly a change in the size of the spectral 
features in the infrared. In this regime, small changes in the 
retrieved planetary radius can compensate for sever al orders 
of ma gnitude of change in the retrieved abundance llGri ffit3 
l2014l) . thus resulting in certain retrieval degeneracies when 
faced with noisy data or data obtained only in the infrared. 

Figure|4]also shows that much of the water features for 
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Figure 7. Transmission spectra without refraction of a Jupiter¬ 
sized exoplanet with a 600 K isothermal Jovian atmosphere for 
various vertical homogeneous mole fractions of water (see inset 
legend). The water-free Rayleigh scattering transmission spectra 
from Fig. [2] are overlaid for comparison (see caption and inset 
legend in Fig. [T] for more details). Refraction will suppress water 
features which lie below the refractive Rayleigh curve associated 
with the host star’s spectral type (coloured triple-dotted dashed 
lines). 
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Figure 8. Transmission spectra with refraction of a Jupiter-sized 
exoplanet, with a 600 K isothermal Jovian atmosphere for vari¬ 
ous vertical homogeneous mole fractions of water (see inset leg¬ 
end), if one could somehow probe the atmosphere to the lower 
refractive boundary, irrespective of the planet-star distance. The 
non-refractive Rayleigh continuum (short-dashed curve) is shown 
for comparison. 


the non-refractive case can lie below the refractive Rayleigh 
continuum, depending on the water abundance and the spec¬ 
tral type of the host star. As these features will be severely 
decreased in the refractive case, Figure [4] illustrates the frac¬ 
tion of water features that can be hidden by refraction for 
various scenarios. In the refractive case, the size of absorp¬ 
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Figure 9. Transmission spectra with refraction of a Jupiter-sized 
exoplanet, with a 600 K isothermal Jovian atmosphere for various 
vertical homogeneous mole fractions of water (see inset legend), 
orbiting a FO star. The non-refractive Rayleigh continuum (short- 
dashed curve) is shown for comparison. 


tion features increases as the critical boundary decreases 
in altitude. The majcimum size occurs if one can somehow 
probe down to the lower refractive boundary, which is shown 
in Figure[Sl Since the effective scaleheight of the atmosphere 
tends to zero toward the lower refractive boundary, refrac¬ 
tion reduces the size of spectral features the closer they are 
to this boundary. This is particularly obvious with the peak 
of the water feature at 2.7 pm for various abundances. With¬ 
out refraction (Figure U)), the size of this feature for a mole 
fraction of 10“^° (green solid line), is similar to the part 
of the feature above the long-dashed solid line for a mole 
fraction of 10“® (blue solid line). However, their sizes in the 
refractive case (Figure^ are completely different. 

What is particularly interesting about Figure [S]is that, 
unlike in the non-refractive case, the size of spectral features 
in the infrared varies dramatically with the abundance of the 
ab sorbing species . Hence, the retrieval degeneracy described 
bv iGriffitH l2014h does not exist for a cloud-free 300 K cool 
atmosphere because the refractive lower boundary provides 
a reference altitude akin to a surface below which spectral 
features cannot exist. Furthermore, the boundary which is 
relevant for the retrieval of atmospheric properties is the 
critical boundary, which raises the refractive continuum and 
decreases the size of spectral features further, and more so 
for hotter stars than for cooler stars. Amongst the spectral 
type that I have considered (see Table [2]), a planet orbit¬ 
ing an FO star has the smallest spectral features. Figure [6] 
shows that around a FO star a 300 K atmosphere with a 
water mole fraction of 10“^'^ does not have any recognisable 
features, and an abundance of 10“® produces only weak fea¬ 
tures, contrary to what is observed for the non-refractive 
case (Figure |4]). 

The breaking of this degeneracy is not due to opacity 
from Rayleigh scattering but rather by two refractive effects 
also responsible for the decrease of the Rayleigh slope: 

• The decrease of the effective scaleheight of the atmo- 
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Figure 10. Transmission spectra without refraction of a Jupiter¬ 
sized exoplanet with a 1200 K isothermal jovian atmosphere for 
various vertical homogeneous mole fractions of water (see inset 
legend). The water-free Rayleigh scattering transmission spectra 
from Fig. [3] are overlaid for comparison (see caption and inset 
legend in Fig. [T] for more details). Refraction will suppress water 
features which lie below the refractive Rayleigh curve associated 
with the host star’s spectral type (coloured triple-dotted dashed 
lines). 
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Figure 11. Transmission spectra with refraction of a Jupiter¬ 
sized exoplanet, with a 1200 K isothermal jovian atmosphere for 
various vertical homogeneous mole fractions of water (see inset 
legend), if one could somehow probe the atmosphere to the lower 
refractive boundary, irrespective of the planet-star distance. The 
non-refractive Rayleigh continuum (short-dashed curve) is shown 
for comparison. 


sphere toward a value of zero as rays approach the lower 
boundary. 

• The deflection, away from the oberver, of rays grazing 
atmospheric regions below the critical altitude. This creates 
a step function with altitude in the perceived illumination of 


Figure 12. Transmission spectra with refraction of a Jupiter¬ 
sized exoplanet, with a 1200 K isothermal jovian atmosphere for 
various vertical homogeneous mole fractions of water (see inset 
legend), orbiting a FO star. The non-refractive Rayleigh contin¬ 
uum (short-dashed curve) is shown for comparison. 


the atmosphere, just as a surface or an optically thick cloud 
does. 

These two effects create a refractive continuum which mani¬ 
fests itself in the infrared, where Rayleigh scatering is weak, 
in spectral regions of low molecular opacity. 

Figures [7}0 and Figures [TCTlll2l display the same infor¬ 
mation as Figures [dHSl but for planetary temperatures of 
600 and 1200 K, respectively. Aside from the scaling of the 
size of the spectral features with temperature, one can ob¬ 
serve only subtle changes with temperature in the overall 
shape of the spectral features in the non-refractive case (see 
Figures 131 [71 and IIOII . In the refractive case, as the plane¬ 
tary temperature increases, differences in the spectral type 
of the host star becomes less important and are almost neg¬ 
ligible at 1200 K (compare Figures ITT] and [T2ll . Also, the size 
of spectral features varies with abundances but these varia¬ 
tions are not as pronounced for the hotter exoplanets as the 
cooler exoplanets (e.g. compare Figures O (El and 1111 when 
one can probe to the lower boundary). This occurs because, 
as the planetary temperature increases, the refractive con¬ 
tinuum is lower relative to the non-refractive H 2 O spectral 
features, and thus has a correspondingly lower impact on 
the shape of these features. 

At temperatures of 1200 K, a water mole fraction of 
10“® or higher produces spectral features which are largely 
unaffected by r e fracti on, a result that seems consistent with 
[Hubbard et al.l (l200lh . However, including refraction allows 
one to constrain the possible range of abundances better. If 
one observes strong features in the data then small abun¬ 
dances are excluded. Conversely, the non-detection of spec¬ 
tral features need no longer be explained only with clouds, 
but can also with low abundances. 

What happens to inferred abundances if an optically 
thick cloud deck is located at a lower density or pressure, 
or higher altitude, than the critical boundary? In this case, 
the cloud deck will raise the continuum and cut down the 
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size of spectral features even further than refraction does, 
so that the spectral features might have the same size as 
if produced by a lower abundance in a clear atmosphere. 
Hence, if one assumes that the atmosphere is cloud-free, 
abundances inferred from transmission spectra will be lower 
limits. 


4 CONCLUSIONS 

I investigated the effects of refraction on transmission spec¬ 
tra of giant exoplanets with a suite of simple isothermal at¬ 
mospheric models for a Jupiter-like exoplanet with a Jovian 
H 2 -He composition, ignoring stellar limb-darkening. The ef¬ 
fects are particularly important for giant exoplanets, more 
so than for terrestrial exoplanets, because they do not pos¬ 
sess a surface which could otherwise limit the density, or 
pressure, to which one can probe their atmosphere. 

Indeed, without refraction, one can probe a pure Jovian 
H 2 -He Rayleigh scattering atmosphere to more than 500 bar 
of pressure at 5 pm, barring other sources of extinction. 
However, refraction creates a thin refractive boundary layer 
at high densities where the perceived density scaleheight 
of the atmo sphere is dramatically smaller th an the actual 
scaleheight (iBetremieux fc Kalteneggerll201 j i. This appar¬ 
ent decrease in the scaleheight reduces the size of spectral 
features and decreases the Rayleigh scattering slope within 
that layer. 

The lower boundary of this refractive boundary layer is 
defined where the refraction-induced curvature of a ray al¬ 
lows it to follow a circular trajectory around the planet. 
Atmospheric regions below this boundary can never be 
probed during an exoplanetary transit. This flattens out the 
Rayleigh continuum in the infrared and effectively mimics 
a surface. The pressure level of this lower refractive bound¬ 
ary increases with temperature and is much lower than in 
the non-refractive case. It varies from a pressure of 5.4 to 
86.3 bar as the atmospheric temperature changes from 300 
to 1200 K in an isothermal Jovian atmosphere. 

The fiat infrared Rayleigh continuum that one ob¬ 
tains with refraction is significantly different from the non- 
refractive case, where the Rayleigh continuum decreases 
with wavelength with a slope characteristic of the scale- 
height of the atmosphere. At shorter wavelengths, the re¬ 
fractive Rayleigh continuum slowly merges with the non- 
refractive one. The wavelength above which refraction 
changes significantly the Rayleigh slope increases with tem¬ 
perature, from about 0.6 to 1 pm between atmospheric tem¬ 
peratures of 300 and 1200 K. Accounting for the angular size 
of the host star seen from the planet, raises the Rayleigh con¬ 
tinuum and changes the Rayleigh slope further, which shifts 
this spectral boundary to lower wavelengths. 

The magnitude of this correction can be quite substan¬ 
tial for an atmospheric temperature of 300 K and for a FO 
spectral type host star, but decreases with decreasing stellar 
temperature and with increasing atmospheric temperature. 
At atmospheric temperatures of 1200 K, the correction for 
a planet around an FO star is very small. At this temper¬ 
ature, the transmission spectrum is almost independent of 
the spectral type of the host star and is mostly determined 
by the location of the lower refractive boundary. 

As I illustrated with water absorption features, the re¬ 


duced effective scaleheight near the lower boundary and the 
step function in the apparent illumination of the atmosphere 
at the critical boundary, also responsible for the reduction of 
the Rayleigh slope, break the retrieval degeneracy between 
abundances and planetary radius. In the non-refractive case, 
when absorption is much stronger than Rayleigh scattering, 
changes of several orders of magnitude in the mole fraction 
of well-mixed species only shifts the spectrum up or down 
in effective radius without changing the size of spectral fea¬ 
tures significantly. However, in the refractive case, the size 
of H 2 O features changes dramatically with abundances for 
a 300 K atmosphere. As the atmospheric temperature in¬ 
creases, the refractive continuum shifts downward relative 
to spectral features, and the variation with abundance of 
the size of spectral features decreases. At atmospheric tem¬ 
peratures of 1200 K, only water mole fraction lower than 
10“® show some significant size variation with abundance. 

The possible presence of an optically thick cloud deck 
complicates retrieval of atmospheric properties further. The 
abundances obtained from assuming a clear atmosphere are 
lower limits. However, the brighter refractive wings just out¬ 
side of transit can be used to distinguish betw een clear and 
cloudy atmospheres (iMisra fc Meadowa 120141 '). Hence, it is 
clear that including all effects of refraction in exoplanet at¬ 
mosphere retrieval algorithms can only help constrain atmo¬ 
spheric properties. 

Although quite important even for hot exoplanets, re¬ 
fractive effects are more pronounced in cooler atmospheres. 
Even though the bulk of exoplanet observations has been 
primarily focused on hot exoplanets because they are easier 
to detect and characterise, long-time monitoring of various 
areas of the sky is just starting to yield longer-period cooler 
exoplanets (ICrossfield et al.ll2015l . Jenkins et aLllioish with 
promises of more to come. Furthermore, with the upcoming 
planned launches of the Transiting Exoplanet Survey Satel¬ 
lite (TESS) and the James Webb Space Telescope (JWST), 
it is only a matter of time before we start obtaining exo¬ 
planet transmission spectra of sufficient quality that signa¬ 
tures of refraction can be exploited to constrain the atmo¬ 
spheric properties of planets of any temperature. 
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